UNCLASSIFIED 







ad 257 797 


iiepAoduced 
ka Ute 


ARMED SERVICES TECHNICAL INFORMATION AGENCY 
ARLINGTON HALL STATION 
ARLINGTON 12, VIRGINIA 




UNCLASSIFIED 













NOTICE: When government or other drawings speci¬ 
fications or other data ere used lor any purpose 
other than in connection with a definitely related 
government procurement operaiicr, the U. S. 
Government thereby incurs no responsibility, nor an> 
obligation vfcatsoever; and the fact that the Govern¬ 
ment may have formulated, furriBhed, or in any way 
supplied the said drawings, specifications, or other 
data is not to be regarded by implication or other¬ 
wise aB in any manner licensing the holder or any 
other person or corporation, or conveying any rights 
or permission to manufacture, use or sell any 
patented invention that may in any way be related 
thereto. 








NASA TN 0-555 


s S'v 

° Qs. 

< f" 

< 

4r» 



i. 


TECHNICAL NOTE 

D-555 


, CD 


XL R OX 




THEORY OP ARTIFICIAL STABILIZATION 
OF MISSILES AND SPACE VEHICLES 
WITH EXPOSITION OF FOUR CONTROL PRINCIPLES 

By R. F. Hoelker 

George C. Marshall Space Plight Center 
Huntsville, Alabama 

a v T \ A 

I A V» * 


W jimiii96ti jH 

liut^ 


1 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON June 1961 








LIST OF SYMBOLS 
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Thru.jt force 
Axial air force 

Air force perpendicular to long axis of 
vehicle 

Control force perpendicular to long axis 
■jf missile., due to motor swiveling or/and 
due to vane ’dcf! potions 

Mass of vehicle 

Moment of inertia of vehicle !.n yaw or 
pitch 

Radius of gyration in yaw or pitch 

Specific, aerodynamic restoring torque 

Specific control torque 

Magnitude of standard velocity of vehicle 

Wind velocity magnitude 

Linear acceleration of center of gravity 
of vehicle perpendicular to standard path 

Velocity of center of gravity of vehicle 
perpendicular to standard path 

Velocity of center of gravity of vehicle 
perpendicular to long vehicle axis 

Local linear acceleration at a. vehicle 
station perpendicular to long vehicle axis 

Angle of attack 

Swivel motor deflection or vane deflection 
Attitude angle 

Wind angle, between flow and standard path 
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ai 

a* 

S* 

Cg 

C P 

C M 

C H 

*i 

Bi 

X 

CD 

c 

Superscripts 

dot 

pi'inio 

Subscripts 

qss 

M 

H 


Attitude displacement gain 

Attitude rate gain 

Attitude acceleration gain 

Angle of attack gain 

Local lateral acceleration gain 

Location of center of gravicy 

Location of aerodynamic center of pressure 

location of accelerometer 

Location of point of interest (as of human 
passenger) 

Coefficients of cnaracteristic equation 
Coefficients of path reaction equation 
Differential operator 
undamped frequency 
Damping ratio 

For differentiation with respect to time 
For differentiation with respect to angle 

For quasi steady state 
For "measured" 

For "at place of interest" as location of 
human passenger 


others 


Explained at place of occurrence 
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THEORY OF ARTIFICIAL STABILIZATION 
OF MISSILES AND SPACE VEHICLES 
WITH EXPOSITION OF FOUR CONTROL PRINCIPLES 

By R. F. Hoeiker 
SUMMARY 

The theory of artificial stabilisation has been developed in a 
sequence of reports, published .it this installation (see References 1, 
2, and 3). The present paper rederives the theory i. a more rigorous 
fashion and extends the application of t‘..e theory beyond the Drift- 
Minimum-Control mode to three further possible modes of control, that 
is the Load-Minimum-Control, the Maximum Comfort Control and the 
Motion Center Control mode. Concern is given in this paper to two 
instrumental arrangements by which the artificial stabilization can 
.materialize. These are the couplings of (a) the space-referenced 
gyro with an angle of attack meter and (b) the space-referenced gyro 
with a missile-mounted acceleromecer. Formulas for both applications 
are derived side by side. The theory is restricted to rigid body 
assumptions. 




INTRODUCTION AND ANALYSIS 


A. DERIVATION OF THE PRINCIPAL EQUATIONS 

Concern is given in this paper to the rigid body motion of a 
space vehicle. Hence, the degrees of freedom introduced by the 
flexibility of the missile or by fluid sloshing are disregarded. 
Further, the reaction of the control members to the coonands of the 
controlling elements is assumed to occur in the idealized form, that 
is, actuator lag, response threshold and inertial effects of the 
swiveling motors are ignored. Further, linearization la applied to 
the representation of the missile motion itself. 

The equations set forth can be understood as representing 
either a vertical missile flight, with the variables measured then 
either in the pitch or yaw plane, or a tilted flight, where the 
variables then reflect the motion perpendicular to the plane of 
flight With the terms as defined in the list of symbols (page ili), 
the equations of motion read: 

Lateral Path Motion: 


z 


X m . N' 

— q> + — a + 



m 
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Angular Motion: 

? + C( a + c 2 ^ * 0 (2) 

Angular Relationship: 

i 

a - a* = 9 - ~ J a w = - ( 3 ) 

For the listing of the control equation, the reader 
will recall from the Introduction that there are several 
modes of materializing the method of artificial stabiliza¬ 
tion, two of which were mentioned. Since it is intended 
here to analyze both, the method based on ang’e of attack 
sensing and that based on missile-fixed acceleA*ometer 
sensing, the analysis is to be spl*c from here on into two 
parallel channels. The formulas to be shown, however, are 
identical with the exception of the terms that concern the 
sensing elements. Therefore, it proves expedient to treat 
the two modes of sensing In the following equations jointly, 
with the understanding that only one of the sensing elements 
will be used in applications. 

The control equation shows the control angle 0, which 
may be a motor swivel angle or the angle of Jet vanes. In 
response to commands from (1st) the attitude 9, as measured 
by reference to space fixed gyros, (2nd) the rate of change 
of cp for damping, (3rd) either the angle of attack a, or the 
local lateral acceleration aw measured by an accelerometer, 
sensitive perpendicular to tne missile's long axis; 

A " *o 9+ *t 9 * £0 « + 8e a M (4) 

I 

V - 

In application, the two 
terms are mutually exclusive. 

The factors a 0 , a t , b Q and gg, are the gain values of the 
control system. Alternately setting later b 0 or g* to zero, 
the effect of the non-zeroed control mode can be studied. 

It is helpful for further reference to write the 
control equation in somewhat modified forms. The acceler¬ 
ation measured by the Instrument can be broken down In two 
parts, that is the linear acceleration, r, of the center of 
gravity of the missile, perpendicular to its axis, and that 
part associated with the angular acceleration about the 
center of gravity. The measured value ag can be written as 

ag - r + ( c m - C 0 ) 9 (5) 
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with Cg and being the locations of the center of gravity 
and the instrument, both measured from the missile cross 
section in which the hinge-point of the control motor is 
located, taken positive in the forward direction. With this 
relationship introduced, the control equation reads: 

Alternate Form of Control Equation: 

0 = a 0 CP + a t $ + b 0 a + Ba p + (Cj| - C a ) <pj (6) 

Sometimes the factors of q> are combined to 

a a 5 - Sft (C M - Co) <p (7) 

which is done to achieve consistency with the gains a 0 and 
a,. The treatment of this term, _a being separate from the 
r - term, is Justified, since by the choice of the measuring 
station cm. the value of a t can be changed independently of 
the value that ga may have. In fact, it is by this choice 
of two parameters that the accelerometer control mode may 
be somewhat superior to the angle-of-attack mode, which 
offers only one parameter (b 0 ) to vary. 

A further substitution may be made, that ties r to the 
angles of the missile motion: 



N* 

— a + 
m 



( 8 ) 


For convenience, the system of equations in the form 
used may be reoeated here. 
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F - X N* R* 

- 0 + — a + — 0 

m mm 


( 9 ) 



( 10 ) 


q> + Cj a + Cg 0 = 0 (ll) 

0 = a o ^ + a 1 ^ + a # »p + b o a+g a (“~a+ 0) (12) 

m m 

By substituting for a and 0 in equations (9) and (11) the 
system may be reduced to the following two equations. 
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Rotary Motion: 


R 1 

qj (1 + a a c a - 8a — ) + ? a i c « + 

ra 


(13) 


N t —| 

+ <p + °2 ( a o + b o) + 8 * (^a ~ " c i ~ 2 | 

■ <t - t> L°' + c * b » + g * (5 *r - ’> i’>] 

and the equation 


5 (1 - * T> + <t - V> 


R' 


.. Ri . R • 

? “i - M a t - + 
m m 


[- 


,, t P - X K* p - 

*■ * K a ° + b 0 - e. — > T + ~' 


x + nH 

m J 
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The left-hand member of Equation (13) is usually con¬ 
sidered describing the dynamical behavibr of the rotary loop 
(i.e., the motion of the vehicle about its center of gravity.) 
We will soon check on the validity of this assumption by 
studying the characteristic equation. 


To arrive at an equation that is informative for the 
reaction of the path motion to the rotary system, we substi¬ 
tute in equation (14) the tera (4 - w)/v by using equation 
(13). This leads to the following equation: 


Path-Reaction: 
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q> 


N' 

R» 

, N' 

R' 


- b 0 — 

- a* (c t — 

-Ci — 
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c i + c t b o + s* (c 8 
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> N' 

C 1 + C 2 b O + (c 8 —“ 

in 



(15) 
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c t 
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N' 

+ c, b 0 + fSm (ca “ 
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Cl 




Ci 
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+ 



+ c a b 0 + & 



which is late? used in the abbreviated fora as 

ZxqpBg + (ji Bj + {pB 0 . 

Before proceeding in the analysis, we list in the following 
the quasi-steady state solutions for a, 0 and f that result, 
if first and higher derivatives of these variables are 
assumed to be zero. A quasi-steady state is approached, if 
the rotary loop is well damped and during a cycle of this 
loon's motion the change of wind and of the lateral linear 
velocity z is negligible. Then e.g. we can directly derive 
from equation (13) the following relationship: 
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^qss 




/ N* 

«. + «■ »0 + *. < c . T - 



(16) 


+ e. 


( a o 


+ b 0 ) + g a (c a 



c i 
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By using equations (10), (11) and (12), we find the following 
quasi-steady states for a and £. 


a qss - 




- a o c* 


(17) 


Cl 


+ c a (a 0 


+ b Q ) + g fl (C 




Pqss 



a 0 ® i 


(18) 


c, + c a 



+ b u ) + g« (c a 



m m 


The developed equations (13)> (15) and (1 6 ) to (13) 
will ho discussed in more detail, after in the next chapter 
the characteristic equation is derived. 


B. THE CHARACTERISTIC EQUATION 

Introducing the operator X and higher powers of it for 
representing the first and higher derivatives of variables, 
the homogeneous part of the system (9) to (12) can be written 
in the following fora: 


z 

1 

c _1 

/3 


- X 

N' 


R' 

P - X 


m 


m 

u 

0 

Cl 


c a 

X* 

1 

v 

1 


0 

- 1 

0 

b o + 2a 

W 

a 

- 1 + g* — 
m 

a o + a t X+ a « X* 

» 


(19) 




From this, the following characteristic equation is developed. 

A 3 X s + A, X* + A. X + A 0 = 0 with 
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The system, in general, has three roots. Though 
principally all three participate in both the path motion 
and the rotary motion, one of the roots can frequently be 
thought of as describing the path motion, while the remain¬ 
ing two aie thought representing the rotary motion. The 
association of motions with roots is approximately valid, 
if the roots are well separated in the complex plane and 
simultaneously the "path-root" is approaching zero. If 
this is the case, then the path root is approximated well 
by the ratio A 0 /A,, where A 0 and A t are the last two terms 
of the characteristic equation. Following this line of 
reasoning, A 0 equal to zero represents a path motion coordi¬ 
nated to a zei'o-root, which is a motion essentially deter¬ 
mined by the initial conditions and constant ( .n velocity) 
with time. This is discussed in more detail in the next 
chapter. 

C. THE DRIFT-MINIMUM-PRINCIPLE 


, By an appropriate choice of the gain values a Q , b 0 
(g* respectively), the zero-order term Ao in the charac¬ 
teristic equation (20) can be made identically zero. As 
we have discussed in the last chapter, this results in one 
root of the system to be zero, which root we associated to 
the lateral path motion z. The physical significance of 
A 0 being zero may more clearly be apprehended, if it ie 
noted that the terra A 0 is recognized also in the "Path- 
Reaction" equation (lf>) representing the coefficient of 
with the exception of the factor v. The structure of (lj) 
is 

z “3iT + Bi<p+Bo<p (21) 

if A 0 - 0, then Bo - 0, which latter condition Is Identical 
to stating that for any steady state condition of q>, that 
is, if $ and £ both are equal to zero,' there is no lateral 
acceleration (z) on the center of gravity of the vehicle. 

Since equation (15) was derived without making any condi¬ 
tion on the wind magnitude, the zero lateral acceleration 
condition la also Independent of the wind magnitude. This 
condition was in the past designated the condition of the 
"Drift-Mlnlmum-Prlnciple." Since in transient motion, 
where $ an<\,$ are of finite values, also transient non-zero 
values for z are to be expected, the claim for zero-drift 
cannot be made* The condition may be written down explicitly: 

Drift-Minimum-Condition 


Lateral linear accelerations from steady-state angles will 
be zero, if the gains a 0 and b 0 (ga respectively) are such* 


that 


c» 

c a 


+ bo + 




c 1 R\ 

Ca tt / 



ao 


F - X 


# 


« 
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( 22 ) 
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The condition requires the two gain values a 0 and bo 
(or a 0 and &) to satisfy only a linear relationship. A 
second condition may be formulated to determine their values. 

The straight line in the (a 0 , b 0 )-plane. that represents 
the linear relationship, divides the (a 0 , b 0 )-plane into two 
areas. Depending on the area in which a combination of ao» b 0 
is situated, the vehicle will fly into the wind or drift with 
the wind. For the first case, the root that is coordinated 
to the path-motion is positive and therefore unstable, for 
the second case,the path root is stable and indicative for 
the rate with which the vehicle's lateral velocity approaches 
the wind velocity. 

The angular attitudes in the Drift-Minimum-Case are true 
steady states a..a follow (16) to (Id) with i being constant 
and small compared with the wind velocity w. The drift- 
minimum case may also be defined as that case where, by the 
control mode, an attitude and throw-angle combination la 
enforced that leads to the cancellation of the sum of all 
force components perpendicular to the nominal flight plane. 

D, THE LOAD-MIMIMUM-PHINCIPLE 

From the equations (17) and (18), it is evident that 
steady-state angles of attack and steady state throw angles 
are reduced by reducing the gain a 0 . An extreme oaae la 
obviously given by having a 0 equal to zero. Since the two 
steady state angles here are zero, thla condition may be 
called the load-minimum-condition. In this situation, the 
missile axis aligns in the flow direction like a weathercock 
stable missile without special attitude control. The atti¬ 
tude angle <p, therefore, can assume large values. A true 
steady-state would be reached only with the flight direction 
being parallel and opposite to the wind direction. Though, 
with respect to the nominal flight plane, the path-motion 
is unstably short time application of this load-minimum- 
control mode for the benefit of load reduction and throw 
angle reduction is not objectionable, provided a temporary 
deviation from the nominal flight plane ia permissible. 

It may be cautioned that the effect of load reduction 
is shown here only for quasi-steady states. No deductions 
can be made from equations (17) and (18) as to the transient 
dynamical behavior of the angles. 

E. THE MAXIMUM COMFOhT PRINCIPLE 


A third principle may briefly be listed, that can be 
effected by means of particular gain value choice. It wil3 
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soon be apparent why this is called the ''maximum comfort 
principle." 

The combination of a linear acceleration ( 2 ) of the 
center of gravity and of an angular acceleration ($) about 
the center of gravity determines for the points along the 
missile axis a linear distribution of the local point 
acceleration, a3 tor instance, measured by local accelerome¬ 
ters, mounted missile-fixed and sensitive perpendicularly to 
the long missile axis. The local normal acceleration ajj 
may be expressed as 


an = t + (Cfj - Cg) q> (23) 


with Cjj - Cq » distance of point of interest (Ch) from 
center of gravity (Cg) positive forward 
of Cg. (The letter "H" taken to indicate 
e. g., the location of a "human passenger.") 

It is conceivable that it is desired for a certain 
point of the vehicle to be free of normal acceleration as 
much as possible. This desire is expressed by the require¬ 
ment that 


a K = ^ T ' "* C (^4) 


Making use of (3) and ( 9 ), equation (24) can be written 
as 


•• 

z 


n v 

r - a 



<P - (Cjj - Cq) <p 


(25) 


which, combined with equation (15) gives the following 
criterion: 
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Maximum-Comfort-Principle 


•• • - _ \ 

m - C.»J - 

' ' *A \J * 


N* R' N' R* 

~ + b 0 — + ag (c a -c, —) 

mm mm 

N' R' 

Ci + c 8 b 0 + g, (c 8 -c t —) 

m m 


N* R* 

_ a t (eg-c, —) 

m m 

mi R* 

ct + ^ b c + g a (eg-c, -) 

m m 


N' R' 

- a 0 (eg- c, -) 

m m 

N' R» 

Ci + Cg b 0 + g, (<*-e, -) 

m m 


A full realization of this criterion v/ould be achieved only 
by having a~ =* ## ai «= 0 and choosing b 0 (or a*, g*) so that 
the factor to y is identically zero. Since cutting of a 0 
and to zero ivill, in general, not be feasible, a partial 
fulfillment of equation [26) only seems to be possible by 
having the factor to $ identically zero. This is expressed 
by the condition: 

[c, (C H - C Q ) - ~ + b 0 [c. (C H - C Q ) - + (27) 

N' R> r -| 

+ (Cg —- - Cj ) gg (Cr - C 0 ) - a # I a 0 
m m [_ _j 
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The true realization of criterion (26) could be of Interest 
e.g., for maintaining a fluid level always perpendicular to 
the long axis of the vehicle. 

F. PRINCIPLE OP MOTION CENTER CONTROL 

As last case, a modification to the Maximum Comfort 
Principle may be mentioned. This is concerned with the 
effort of controlling the point of minimum motion, which 
is the point on the long axis of vehicle, at which the rota¬ 
tion about the center of gravity and the linear motion of 
the center of gravity perpendicular to the nominal flight 
plane cancel each other. In contrast to the former case, 
here the acceleration 2, perpendicular to the nominal path, 
instead of r, which is normal to the missile axxs, is to be 
considered. The acceleration of * ’-'ehicle point (Ch - Co), 
perpendicular to the nominal path, is expressed by 

5 h •• 5 + (Ch - Co) q>, which for z# - 0, leads to 
z « - (Ch - Co) f. (28) 


Feeding this requirement into the path-reaction equation 
(15) results in the following criterion: 

Motion Center Control 




N» R* N' R» H 

— + b 0 — + a a (c t -c x —) 

mm m m | 


N» R' 

c i + c a *o + £e Us- c i —> 

m m J 

N« R' 

- a, (c # — - Cj —) 
m m 

N» R' 

Ci + c* b Q + gfc (c g -c, —) 

m m 


N« R* F - “X r~ 

- *o (c» - - c i ' ) + - c x + c, 

mm m 


M» R* 

b o + Be ( c * - - c* ” 
m m 
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It is interesting tonote by comparison of (29) with (26) 
that thC factors of qj and <p are identical in the two equa¬ 
tions, that, however, the factors of <p in ( 26 ) and ( 29 ) 
reflect exactly the two earlier discussed principles of 
load minimum and drift minimum control. 


me pnysical significance of the Notion Center Control 
Mode is that at the location of the Motion Center, e. g., a 
fluid layer remains parallel to the horizontal (in vertical 
flight) and the layers above the motion center experience 
an acceleration opposite in sign to those below the motion 
center. (This assumes no interaction between layers.) 

If we ignore the term with <p, the inter-rel'tionship 
between the four discussed principles of control nay then 
be summarized in the following way; 


1. The minimum load principle calls for a 0 - 0. 


2. The minimum drift principle requires that a<> and 
bo (or g a ) satisfy equation ( 22 ) which reads: 


Cl 

c* N' Cj R* N' - 5 - R» 

- + b 0 + g* (-) - a 0 -*— 

c a a c, m P - X 


( 22 ) 


3. The maximum comfort principle (or vehicle related 
control of fluid levels) calls for a combined satisfaction 
of ( 1 st) the load minimum principle and ( 2 nd) the equation 
( 27 ) which reads: 



NH 


(C H - C Q ) - —-J + b 0 (Ch - Cq) - 



+ 



4. The motion center control requires satisfying also 
the Just listed equation ( 27 ), but now in combination with 
the drift minimum principle, (equation 22 ). 

Q. LIMITATIONS IN SATISFYING SIMULTANEOUSLY SEVERAL 

-raiiBIHSS- 


The possibility of satisfying several requirements 
simultaneously depends on the number of open parameters that 
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are available. The variable ai Is determining primarily the 
degree of damping of the rotary motion, as shown in the 
<p-term of equation (13), though the rooting of equation (20) 
shows that the damping Is also slightly dependent on other 
parameters The damping characteristic is usually allowed to 
vary in large limits (as from 4o£ to 140# of critical damping), 
but it is not allowed to be zero. This excludes the term of 
<p in (15) to vanish identically. 

A demand on ao and bo (or g«) is imposed from consider¬ 
ations of the rigid body frequency. An approximate expres¬ 
sion for the frequency of the rotary motion is derived from 
the coefficient of $ in equation ( 13 ) with 


N' R' 

Ci ■» c a (a 0 + b 0 ) + gs (c a - - Ci —) 

m m 


R' 

1 + a a c a - g* — 

m 


(30) 


with a: being the undamped frequency (in rad per sec). 

It Is obvious that two requirements can be met, since 
two gains are at our disposition, a 0 and bo (g«, respectively). 
Thus, besides the frequency, we can e. g., choose to go to 
drift minimum or load minimum or any path reaction mode 
between these two. 

Further, in the particular v.ase of utilizing an acceler¬ 
ometer, the location of this Instrument along the long vehicle 
axis can be chosen arbitrarily (In bounds) which provides us 
a third variable. This location is reflected in the magni¬ 
tude of the gain a s . By virtue of this variable, the 
accelerometer-based vehicle control arrangement is superior 
to the angle of attack-based arrangement. (It may be 
recalled nere that concern is given here only to the rigid 
body mode.) 

H. COMPILATION OF FORMULAS IN REDUCED FORM 


The formulas listed in the previous chapters adhere 
to the symbols that were introduced in the original force 
and moment equations. Ey taking account of the relationships 
that exist between forces and moments, the formulas will be 
simplified, which also facilitates better understanding for 
the relative significance of the physical parameters. 






2Q 


15 


Distances along the long vehicle axis are referred to 
the cross-plane of the motor swivel points as origin, this 
location (if needed) being marked by C(j. 

Cp - location of "control center", that is, hinge 
point of motor or of control vanes. 

- location of center of gravity (for pitch or yaw 
motion). 

Cp ■ location of aerodynamic center of pressure of 
vehicle 

Cji «■ location of accelerometer, measurii^ perpendicularly 
to the long missile axis 

Cu » location of point of interest, as that of the human 
pilot 

The direction forward of tne control center, Cc $ ia 
assigned the positive sign. Since Cc is assigned the origin, 
the relations hold: 


C G ” C C ” C Q 
Cp - Cc * Cp, etc. 

Besides these signed linear dimensions, the only other linear 
quantity that Is Introduced is the radius of gyration r, which 
18 positive and defined by 

e * m r* 

where e is the moment of inertia In pitch or yaw and a ie the 
mass of the vehicle. 

Then the aerodynamic restoring torque is given by 
dN 

— a (C 0 - C P ) - N» a (Cq - Cp) 
da 

from which follows for c t : 

N' (Cq - C P ) 


The control torque (from motor or vane-deflections) is 
given by dR 

~ A C G = £ C q 

d0 
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from which follows: 



m r a 


Further, the ratio c t /c t is expressed by: 


Ci N‘ (C Q - C P ) 

_ m ■ i . . 

C» R' C 0 


The following terms that frequently occur are listed for 
convenience, using the just derived relationships: 

N* rl' r “1 N' R» 

c a N' - c, R 1 <= - Cq - (Cq - Cp) • - Cp 

mr* L m r* 


Cl 

c a 


R* 


- IV 


Cq - Cp 
” c 0 


N* (1 


C 1 

c C 


) 


c i Cp 

N'-R' - N' — 

c » Cq 


Using these relationships, we reforaulate the most Important 
of the formulas derived in the fonaer chapters. 





IT 


juatior 


9 (1 - a a 


- Se —) + 9 a 


. «' Cq 


(13-A) 


m r* 


P* (Co - Cp) , % N R' C Q g, N» R» m ~[ 

+ <P - + ( a o + b o) -* +-- Cp 

L m r* n r* m m r* 

Z w fN» (Cq - C P ) R' C Q g, N' rt' “| 

■ (-) -;- + b o —r + - a c p 

v vjmr* m r* m m r 8 "J 

Undamped frequency of Rotary Motion, if sufficiently 
separated from the path motion: 

(a) for angle of attack control: 

L (c 0 . Cp) + R' C 0 (a 0 + b 0 N (30-A) 

m r* L J 


(b) for accelerometer control: , N i 

N» (C 0 - Cp) + R» (a 0 C a + g* — Cp) 


1 + gt 


Z (*jJ° . A 

» \Sa^ / 


(30-B) 


Using equation (7), the last formula can be written by the 
use of Cm which is the location of the measuring accelerome¬ 
ter, instead of using a t : 


N' (C Q - C P ) + R' (a 0 C Q + g, - Cp) 

1 m 


m r* R 1 /i 

1 + *t( 


Cq (Cm - Cq) 


- 1 ) 


(30-C) 



18 


Ratio £ of damping to critical damping: 

(a) for angle of attack control: 

dj (31-A) 

m r* 2 
*•« « 

c -- 

2 (o 
or 

C = a * R ‘ Cq (31-B) 

2^ m r*[jj' (Co - C P ) + (a'o + b 0 ) R' CqJ 

(b) for accelerometer control: 


C 


x= 


a t R 


• r. s. 


m 


Lr« 


Cq (Cjh - Cq) 



(31-C) 


1/ F N* T 

y m r* N* (Cq - Cp) + R‘ (a 0 Cq + g* — C p ) 

L m _ 


Path Reaction (equation 15) 

(a) for angle of attack control: 


z = <p r* 


N’ + b 0 R* 


R' 

&x - N' Cp 


■+<p 


N» (C p -Cq) - b 0 R* Cq N» (C p -Cq) - b 0 R* Cq 


R' 

ao - 5 - N« C P 

N* (Cp - C Q ) - b 0 R' Cq 



m 


(15-A) 
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(b) for accelerometer control: 

R* 

r* + a„ — C r 


w p 
m 


z = 


<P 


R' 

Cp - Cq •• g t - Cp 

m 


R‘ 

»i C P 
m 


+ ¥ 


R' 

Cp - Cq — gt Co 
m 


(15-B) 


* 


* 


+ 


R* 

a 0 — Cp 
m 


9 


R» 

C P - Cq - g, - Cp 

a 



m 


Again in the last form, at can be replaced by 
Sm (Cm - Cq), to show directly the influence of the loca¬ 
tion of the measuring instrument: 


” H 

2 » <p 


r* + Et 


R« 

~ Cp (C M - Cq) 


R' 

Cp - Cq - gt — Cp 
m 


R' 

a,— C p 
a 

+ - + 

R' 

Cp - Cq - gt — Cp 

m 


+ 


¥ 


R» 

a 0 — Cp 


R' 

Cp - Cq - gt —— Cp 
m 



(15-C) 


> 
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Quasi-Steady States (equations 16 thru 18) 

(a) for angle of attack control: 

- N' (Cg - Cp) - R' C Q b 0 w - z 
<Pq3s N' (Cq - Cp) + R' c Q (a 0 + b7) V 

a 0 R* Cq W - z 

C * qSS N' (Cq - C P ) + R' Cq (a 0 + b 0 ) V 

- a 0 N' (Cq - Cp) W - z 

Anss = —.. 11 —. ... . 1 

N* (Cq - C P ) + R« Cq (a 0 + b 0 ) v 

(b) for accelerometer control: 


(l6-A) 


(17-A) 


(18-A) 


Oqss - 


“ R' “ 

N* Cp - Cq - gt — Cp 
- m — 


w - z 


R' 

V 

a 0 R* Cq - N* (Cp - Cq - g. 

—- 

Cp) 


m 


ao R 1 Cq 


• 

w - z 


R» 

V 

a 0 R' Cq - »?' (C P - Cq - g* 

— 

Cp) 


m 


- a 0 N« (Cq - Cp) 


w - z 


(l6-B) 


(17-B) 


(18-B) 


R» Cq - N 1 (Cp - Cq- g, - C p ) 


I 
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The Drift-Mnimum-Condition (equation 22) 

(a) for angle of attack control; 

Cp Cq Cq - Cp 

- _ b 0 - ______ 

F - X N‘ R' 


(22-A) 


(b) for accelerometer control; 


Cp ^ Cp 

—— “ Kq —. 

' F - X “ m 


C„ - Cl 

u l 


(22-B) 


The supplementary condition for both the Maximum Comfort 
Principle and Motion Center Control Mode.(equation 2?) 


(a) for angle of attack control: 

N* < C Q “ C P) ( C H - Cq) - r « 
«' Cq (C H - C 0 ) - r» 


b 0 - - 


(27-A) 


(b) for accelerometer conti'ol: 

m (C 0 - C P ) (C H - C G ) - v* 


(Ch - C Q ) - a a 


C P 


(27-B) 


When substituting for a a according to equation (7): 
m (Cq - Cp) (Cji - Cq) - r* 


8 « - + 


R' 


(Cm - c H ) Cp 


(,27-C) 
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